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Abstract 

MACRO was a multi-purpose experiment that took data from 1989 to 2000, at the underground Laboratory of 
Gran Sasso (Italy). MACRO gave important results/limits concerning: (i) the oscillation of atmospheric neutrinos, 
also in the non-conventional scenario of violations of Lorentz invariance, (ii) the searches for exotic particles 
(supermassive GUT magnetic monopoles, nuclearites, WIMPs), (iii) muon physics and astrophysics. A summary of 
the MACRO results will be presented and discussed, focusing the attention on the exotica searches. 



1 Introduction 

MACRO was a large area multipurpose underground detector [1] designed to search for rare events and rare 
phenomena in the cosmic radiation. The detector was located in Hall B of the underground Gran Sasso 
Laboratory (Italy). It was optimised to search for the supermassive magnetic monopoles [2, 3] predicted 
by Grand Unified Theories (GUT). The experiment obtained important results on atmospheric neutrino 
oscillations [4, 5, 6, 7, 8] and performed neutrino astronomy studies [9], indirect searches for WIMPs [10], 
search for low energy stellar gravitational collapse neutrinos [11], studies of the high energy underground 
muon flux (which is an indirect tool to study the primary cosmic ray composition and high energy hadronic 
interactions [12]), searches for fractionally charged particles (LIPs) [13] and other rare particles that may 
exist in the cosmic radiation. 

The detector started data taking in 1989 and it was running until December 2000. The apparatus had 
global dimensions of 76.6 x 12 x 9.3 m 3 and was composed of three sub-detectors: liquid scintillation counters, 
limited streamer tubes and nuclear track detectors. Each one of them could be used in "stand-alone" and 
in "combined" mode. It may be worth to stress that all the physics and astrophysics items listed in the 
1984 proposal were covered and good results were obtained on each of them. 



2 Atmospheric neutrino oscillations 

MACRO detected ^-induced muon events in 4 different topologies. 

The upthroughgoing muons come from interactions in the rock below the detector; the v^'s have a 
median energy of ~ 50 GeV. 

Fig. 1 shows the zenith distribution of the measured 902 upthroughgoing muons (black circles) compared 
with two MonteCarlo (MC) predictions: the Bartol96 [14] flux with and without oscillations (the dashed 
and solid lines, respectively) and the Honda2001 flux [15]. The FLUKA MC predictions [16] agree perfectly 
with the Honda2001. 

For a subsample of ~ 300 upthroughgoing events, we estimated the muon energy through Multiple 
Coulomb Scattering in the rock absorbers in the lower apparatus [7]. The evaluated resolution on E v is 
~ 100%. The parent neutrino path length is L ~ 2i?gcosO, where Re is the Earth radius. Fig. 2 shows 
the ratio data/MC as a function of the estimated L/E v for the upthroughgoing muon sample. The black 
circles are data/Bartol96 MC (assuming no oscillations); the solid line is the oscillated MC prediction for 
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Figure 1: Comparison of the MACRO 
upward-throughgoing muons (black circles) 
with the predictions of the Bartol96 and 
of the Honda2001 MC oscillated and non 
oscillated fluxes (oscillation parameters 
Am 2 = 2.3 • lfT 3 eV 2 and sin 2 29 m = 1). 



Figure 2: Ratio Data/MCno osc as a function of the 
estimated L/E v for the upthroughgoing muon sample 
(black points). The solid line is the MC expectation 
assuming Am 2 = 2.3 ■ KT 3 eV 2 and sin 2 26 m = 1. The 
last point (black square) is obtained from the IU sample. 



Am 2 — 2.3 • 1CP 3 eV 2 and sin 2 26 m = 1. The shaded region represents the simulation uncertainties. The 
last point (black square) is obtained from semicontained upward going muons. 

The low energy events (IU, ID + UGS [5]) are produced by parent Vn interacting inside the lower 
detector, or by upgoing muons stopping in the detector. The median energy of the parent neutrino is 
~ 3 — 4 GeV for all topologies. In both cases, the zenith distributions are in agreement with the oscillation 
prediction with the optimised parameters [8]. 

In order to reduce the effects of systematic uncertainties in the MC absolute fluxes we used the following 
three independent ratios [8]: 

(i) High Energy data: vertical/horizontal ratio, R\ = N vert /Nh or 

(ii) High Energy data: low energy/high energy ratio, R 2 = N low /N high 

(iii) Low Energy data: i? 3 = (Data/MC)w /(Data/MC)io+UGS 

Combining the three independent results, the no oscillation hypothesis is ruled out at the ~ 5c level (6a if 
the absolute values compared to the Bartol96 flux are used) [17]. 

To evaluate the hypothesis of oscillation for different values of Am 2 and sin 2 28 m , the Feldman-Cousins 
[18] procedure was used and the corresponding 90% C.L. region for the < — ► v T oscillation is given in ref 
[8]. The best fit is reached at Am 2 = 2.3 ■ 1(T 3 eV 2 and sin 2 29 m = 1. 

2.1 Search for exotic contributions to atmospheric neutrino oscillations 

MACRO searched for "exotic" contributions to standard mass-induced atmospheric neutrino oscillations, 
arising from a possible violation of Lorentz invariance (VLI) , using two different and complementary 
analyses. The first approach uses the low energy (E v < 28 GeV) and the high energy (E v > 142 GeV) 
samples. The mass neutrino oscillation parameters have the values given in Sect. 2 and we mapped the 
evolution of the \ 2 estimator in the plane of the VLI parameters Av and sin 2 26 v . No \ 2 improvement 
was found, so we applied the Feldman-Cousins [18] method to determine 90% C.L. limits on the parameter: 
\Av\ <3-10- 25 [19]. 

The second approach exploits a data subsample characterised by intermediate neutrino energies. It is 
based on the maximum likelihood technique and considers the mass neutrino oscillation parameters varying 
in the 90% C.L. border [8]. The obtained 90% C.L. limit on the Av parameter is also around 10~ 25 [20] . 
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3 Neutrinos from astrophysical sources 



3.1 Search for astrophysical HE muon neutrinos 

High energy v^s are expected to come from several galactic and extra-galactic sources. An excess of events 
was searched for around the positions of known sources in 3° (half width) angular bins. The 90% C.L. 
upper limits on the muon fluxes from specific celestial sources were in the range 10 -15 -j- 10~ 14 cm~ 2 s" 1 
[21]. A search for time coincidences of the upgoing muons with 7-ray bursts was also made. No statistically 
significant time correlation was found [9] . 

A different analysis was made for the search for a diffuse astrophysical neutrino flux, using a dedicated 
method to select higher energy upthroughgoing muons. The flux upper limit was set at the level of 1.5 • 10~ 14 
cm" 2 s- 1 [22]. 

3.2 Indirect searches for WIMPs 

Weakly Interacting Massive Particles (WIMPs) could be part of the galactic dark matter; they could be 
intercepted by celestial bodies, slowed down and trapped in their centres, where WIMPs and anti- WIMPs 
could annihilate and yield neutrinos of GeV or TeV energy, in small angular windows from their centres. 
One WIMP candidate is the lowest mass neutralino. 

To look for a WIMP signal, we searched for upthroughgoing muons from the Earth centre, using 10° -j-15° 
cones around the Nadir; the 90% C.L. muon flux limits are 0.8 4- 1.4 • 10~ 14 cm~ 2 s _1 [10]. These limits, 
when compared with the predictions of a supersymmetric model, eliminate a sizable range of parameters 
used in the model. 

A similar procedure was used to search for from the Sun: the muon upper limits are at the level of 
about 1.5 2 • 10~ 14 cm" 2 s" 1 [10]. 

3.3 Neutrinos from stellar gravitational collapses 

A stellar gravitational collapse of the core of a massive star is expected to produce a large burst of all 
types of neutrinos and antineutrinos with energies of 5 60 MeV and with a duration of ~ 10 s. No stellar 
gravitational collapses in our Galaxy were observed from 1989 to 2000 [11]. 

4 Search for exotic particles 

4.1 Search for GUT magnetic monopoles (MMs) 

Supermassive magnetic monopoles predicted by Grand Unified Theories (GUT) of the electroweak and 
strong interactions should have masses tom ~ 10 17 GeV. 

MACRO was optimised to search for an isotropic flux of GUT MMs in the cosmic radiation. The three 
sub-detectors had sensitivities in different (3 regions, covering the velocity range 4 • 10~ 5 < (3 < 1. They 
allowed multiple signatures of the same rare event candidate. No candidates were found by any of the three 
subdetcctors. Fig. 3 shows the 90% C.L. flux upper limits for g = go poles (one unit of Dirac magnetic 
charge) plotted versus (3 [2] together with direct limits set by other experiments [23]. The MACRO MM 
direct limits are by far the best existing over a very wide range of (3. 

The interaction of the GUT monopole core with a nucleon can lead to a reaction in which the nucleon 
decays, M + p — > M + e + + tt°. MACRO dedicated an analysis procedure to detect nucleon decays induced 
by the passage of a GUT MM in the streamer tube system (a fast e + track from a slow ([3 ~ 10~ 3 ) MM 
track). The 90% C.L. flux upper limits established by this search are at the level of ~ 3 • 10~ 16 cm~ 2 s _1 
sr _1 for 10~ 4 < (3 < 0.5 • 10~ 2 ; they are valid for catalysis cross sections 5 • 10 2 < <r cat < 10 3 mb [3]. 

4.2 Search for nuclearites, Q-balls and LIPs 

Strangclcts should consist of aggregates of u, d and s quarks in almost equal proportion [24] and would have 
typical galactic velocities (3 ~ 10~ 3 . The MACRO 90% C.L. upper limits for an isotropic flux of nuclearites 
with 10~ 5 < (3 < 1 was at the level of 1.5 • 10~ 16 cm~ 2 s^ 1 sr^ 1 [25]. 
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Figure 3: 90% C.L. upper limit obtained by MACRO 
for an isotropic flux of GUT MMs with g = 
go compared with direct limits given by other 
experiments. 



Figure 4: 90% C.L. upper limits versus mass for 
downgoing nuclearites with (3 = 2 • 10 -3 at ground 
level. The MACRO limit for nuclcarite masses 
larger than 5 • 10 22 GeV/c 2 has been extended and 
corresponds to an isotropic flux. 



MACRO searched also for charged Q-balls (aggregates of squarks, sleptons and Higgs fields) [26], giving 
an upper limit of ~ 10~ 16 cm~ 2 s _1 sr^ 1 [27]. 

Fractionally charged particles could be expected in Grand Unified Theories as deconfined quarks; the 
expected charges range from Q — e/5 to Q = 2/3e. LIPs should release a fraction (Q/e) 2 of the energy 
deposited by a muon traversing a medium. The 90% C.L. flux upper limits for LIPs with charges e/3, 2/3e 
and e/5 are at the level of 10~ 15 cm~ 2 s _1 sr _1 [13]. 



5 Conclusions 

— Standard atmospheric oscillations: no-oscillation hypothesis ruled out at 5 4- 6a. 

Best fit parameters: Am 2 = 2.3 • 10~ 3 eV 2 and sin 2 26 m = 1. 

— VLI: A„ upper limits of the order of 10 -25 . 

— MM search: upper flux limit of 1.4 ■ 10~ 16 cm -2 s -1 sr -1 for 4 • 10 -5 < < 1. 

— Nuclearite search: upper flux limit of 10~ 16 cm~ 2 s _1 sr -1 for f3 ~ 10 -3 . 

— Charged Q-balls search: upper flux limit of ~ 10" 16 cm~ 2 s _1 sr -1 . 

— WIMP search: upper flux limit of ~ 10~ 14 cm~ 2 s _1 sr -1 . 

— LIP search: upper flux limit of 6.1 • 10~ 16 cm~ 2 s _1 sr -1 . 
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